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ABSTRACT 

We describe a mechanism by which supermassive black holes can form directly in the 
nuclei of protogalaxies, without the need for "seed" black holes left over from early 
star formation. Self-gravitating gas in dark matter halos can lose angular momentum 
rapidly via runaway, global dynamical instabilities, the so-called "bars within bars" 
mechanism. This leads to the rapid buildup of a dense, self-gravitating core supported 
by gas pressure — surrounded by a radiation pressure-dominated envelope — which 
gradually contracts and is compressed further by subsequent infall. We show that 
these conditions lead to such high temperatures in the central region that the gas 
cools catastrophically by thermal neutrino emission, leading to the formation and 
rapid growth of a central black hole. 

We estimate the initial mass and growth rate of the black hole for typical condi- 
tions in metal-free halos with T v i r ~ 10 4 K, which are the most likely to be susceptible 
to runaway infall. The initial black hole should have a mass of < 20M Q , but in prin- 
ciple could grow at a super-Eddington rate until it reaches ~ 10 4 — 10 6 M Q . Rapid 
growth may be limited by feedback from the accretion process and/or disruption of the 
mass supply by star formation or halo mergers. Even if super-Eddington growth stops 
at ~ 10 3 — 1O 4 M , this process would give black holes ample time to attain quasar-size 
masses by a redshift of 6, and could also provide the seeds for all supermassive black 
holes seen in the present universe. 

Key words: black hole physics — cosmology: theory — galaxies: formation — ac- 
cretion, accretion discs — instabilities — hydrodynamics 



1 INTRODUCTION 

Several scenarios have been presented for the formation and 
growth of supermassive black holes (SMBHs) in the nuclei 
of galaxies. One possible route traces the black-hole progen- 
itors back to the first generation of stars. The first stars 
formed out of metal-free gas, with the lack of an efficient 
cooling mechanism possibly leading to a very top-heavy ini- 
tial stellar mass function (Larson f998; Carr, Bond & Ar- 
nett 1984). Numerical simulations of the fragmentation of 
primordial clouds in standard cold dark matter theories sug- 
gest that Pop III stars were indeed very massive (Bromm, 
Coppi & Larson 1999, 2002; Abel, Bryan & Norman 2000), 
and would have left behind black hole "seeds" of anywhere 
from 10's to several hundred solar masses. The main features 
of a plausible scenario for the hierarchical assembly, growth, 
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and dynamics of massive black holes from such seeds have 
been discussed most recently by Volonteri, Haardt & Madau 
(2003), Volonteri et al. (2005), and Volonteri & Rees (2005). 

Another family of models for massive black hole forma- 
tion is based on the collapse of supermassive objects formed 
directly out of dense gas (Haehnelt & Rees 1993; Umemura, 
Loeb & Turner 1993; Loeb & Rasio 1994; Eisenstein & Loeb 
1995; Bromm & Loeb 2003; Koushiappas, Bullock & Dekel 
2004) . The main challenge for these models is the disposal of 
angular momentum. Eisenstein & Loeb (1995) and Koushi- 
appas et al. (2004) investigated the formation of black holes 
from low angular momentum material, either in halos with 
extremely low angular momentum (Eisenstein & Loeb 1995), 
or by considering only the low angular momentum tail of ma- 
terial in halos with efficient gas cooling. But even in these 
models, as in all the others, substantial angular momentum 
transport is required in order for the gas to form a cen- 
tral massive object, which ultimately collapses as a result 
of the post-Newtonian gravitational instability. Various an- 
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gular momentum transport mechanisms have been invoked, 
including radiation drag against the cosmic microwave back- 
ground (at very high redshifts: Umemura et al. 1993), vis- 
cosity driven by magnetic fields or turbulence, Rossby waves 
(Colgate et al. 2003) and self-gravitational instabilities. 

The scenario we investigate here is related to the 
second family of models, and focuses on the outcome of 
global dynamical instabilities driven by self-gravity, the so- 
called "bars within bars" mechanism (Shlosman, Frank & 
Begelman 1989; Shlosman, Begelman & Frank 1990). Self- 
gravitating gas clouds become bar-unstable when the level 
of rotational support surpasses a certain threshold. A bar 
can transport angular momentum outward on a dynamical 
timescale via gravitational and hydrodynamical torques, al- 
lowing the radius to shrink. Provided that the gas is able 
to cool, this shrinkage leads to even greater instability, on 
shorter timescales, and the process cascades. This mech- 
anism is a very attractive candidate for collecting gas in 
the centers of halos, because it works on a dynamical time 
and can operate over many decades of radius. In contrast to 
the formation of a supermassive "star", with high entropy 
throughout, we show that the "bars within bars" mechanism 
produces a "quasistar" with a very low specific entropy near 
the center. As a result, the initial core collapse leading to 
black hole formation involves only a few solar masses, rather 
than the thousands of solar masses usually associated with 
direct collapse models. Despite this modest beginning, ac- 
cretion from the envelope surrounding the collapsed core 
can build up a substantial black hole mass very rapidly — 
possibly at a highly super-Eddington rate. 

The plan of the paper is as follows. In § 2, we discuss the 
criterion for global gravitational instability and apply it to 
the gas in dark matter halos with a realistic distribution of 
angular momentum parameters. If more than a few tenths 
of the baryonic matter fall toward the center of the halo, 
then gravitational instability should be very common. But 
even an infall fraction of ~ 10% can lead to an interesting 
number of unstable halos. In § 3 we specialize to halos with 
virial temperatures T v ; r > 10 4 K and metal-free gas. The 
bars within bars scenario makes specific predictions about 
the radial distribution of infalling gas and the associated cir- 
cular velocity, which goes from constant in the outer parts of 
the inflow to quasi-Keplerian close-in. We first discuss the 
infall processes neglecting star formation, and then show 
how the process is modified (but not necessarily halted) if a 
fraction of the infalling gas at each radius forms stars. The 
gravitational binding energy liberated by infalling gas in- 
creases steadily with decreasing radius, until the luminosity 
exceeds the Eddington limit, the infall stalls, and a radia- 
tion pressure-supported "quasistar" forms (§ 4). The radius 
of the quasistar is a few astronomical units, a scale that does 
not change even as the quasistar grows in mass. 

We show that the quasistar has a positive specific en- 
tropy gradient, and that gas pressure remains important in 
the quasistar's core. The temperature of this core steadily 
increases as matter piles on to the quasistar (§ 5), eventually 
approaching 10 9 K, at which point it undergoes catastrophic 
cooling and collapse by thermal neutrino emission (§6). We 
argue that this leads to the formation of a black hole of 
~ 10 — 2OM , which may then grow at rate that greatly 
exceeds the Eddington limit (§7). This rapid growth could 
produce a black hole of several million solar masses, although 



feedback and depletion of the mass supply could quench the 
growth rate at an earlier stage. We discuss the co-evolution 
of the black holes and their hosts and the global impact of 
the black hole population in § 8. We conclude by discussing 
the implications of this model for the interpretation of high- 
z quasars, the statistics of black hole masses in the local 
universe, and its relevance to other astrophysical situations 
where black holes could grow at a very rapid rate. 

Unless otherwise stated, all results shown below refer to 
the currently favored ACDM world model with Om = 0.3, 
= 0.7, h = 0.7, Q b = 0.045, cr 8 = 0.93, and n = 1. 



2 CONDITIONS FOR RUNAWAY COLLAPSE 

We focus here mainly on the dynamical stability of the gas 
in halos with virial temperatures T v i r > 10 4 K. Runaway col- 
lapse could also occur in smaller halos, provided that molec- 
ular hydrogen cooling is efficient and gas can cool well below 
the virial temperature. In the absence of molecular hydro- 
gen, gas in halos with T v i r < 10 4 K would tend to remain less 
dense than the dark matter; tidal forces would then prevent 
widespread collapse and fragmentation at this stage. Since 
cooling and collapse of the gas is more likely in large halos, 
and the masses involved are larger, we henceforth refer to 
halos with virial temperatures T vlI > 10 4 K , unless other- 
wise stated. We stress nevertheless that runaway collapse is 
not completely ruled out in smaller systems at early times, 
well before the first generation of stars created a photodis- 
sociating background. 

Bromm & Loeb (2003) show that if molecular hydro- 
gen formation is suppressed in halos with T v i r > 10 4 K, the 
gas tends to condense into massive clumps in the center. 
The gaseous component of these halos can cool even in the 
absence of H2 via neutral hydrogen atomic lines to ~ 8000 
K, and contract nearly isothermally (Oh & Haiman 2002). 
These massive clumps do not fragment as long as molecu- 
lar hydrogen remains unimportant. One way to hinder the 
formation of molecular hydrogen is the presence of a disso- 
ciating background (Haiman, Abel & Rees 2000; Bromm & 
Loeb 2003; Oh & Haiman 2002; but see Machacek, Bryan 
& Abel 2001). It is therefore possible that the formation of 
seed black holes in massive halos follows an earlier epoch 
of star formation. We consider redshifts high enough that 
a large fraction of gas is still unenriched by metals, or very 
lightly polluted, so that metal line cooling is still unimpor- 
tant (Santoro & Shull 2006). 

We assume that the baryons preserve their specific an- 
gular momentum during collapse (Mo, Mao & White 1998), 
and settle into a rotationally supported disc at the center 
of the halo (Mo et al. 1998; Oh & Haiman 2002). Flattened 
systems can be subject to dynamical and secular instabil- 
ities, even when embedded in external halos (Fall & Efs- 
tathiou 1980). Several instability criteria have been investi- 
gated (e.g., Ostriker & Peebles 1973; Efstathiou, Lake & Ne- 
groponte 1982; Christoudoulou, Shlosman & Tohline 1995), 
which determine the maximum rotational energy (or angu- 
lar momentum) that a system can possess and still be stable 
against bar- like instabilities. Christoudoulou et al. (1995) 
propose a simple, but robust, criterion for stability which 
can be expressed as 
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where T is the rotational kinetic energy, W is the gravita- 
tional potential energy, and / is a parameter dependent on 
the geometry of the system, with / = 1 for discs. 

Numerical simulations have not yet clarified the detailed 
dynamics of gaseous collapse in halos, and we explore here 
three different models for self-gravitating gas discs. We as- 
sume that the disc has either constant circular velocity (Mes- 
tel discs: Mestel 1963) or constant angular velocity (rigid 
body rotation), or that the gas settles down into a clas- 
sical exponential disc. We embed the gaseous discs into a 
halo of mass Mh, virial radius R v i T , and virial temperature 
T v i r , described by a Navarro, Frenk & White (1997, here- 
after NFW) dark matter density profile, with a spin param- 



eter A sp in (= JhEj/ 2 /GM^'^ , where Jh is the total angular 
momentum and Eh is the binding energy). We recall that, 
within the spherical collapse model, the mass of a halo, at 
a given redshift of formation, scales with the virial temper- 
ature as M h ~ 10 4 A~ i 1 r /2 T~^ /2 M e , where A vir is the virial 
density in units of the critical density. 

We determine the characteristics of the gaseous discs 
via a procedure similar to that of Mo et al. (1998). We then 
apply the Christoudoulou et al. (1995) criterion in order to 
determine the stability of the modeled systems. Stability 
depends on two parameters, the halo spin parameter A sp m, 
and the fraction of baryonic matter that ends up in the disc, 
fd = (£lM/Qb)(M<ii B c/Mh)- The results can be understood 
qualitatively by approximating the disc kinetic energy as 
Tdisc ~ 0.5 Mdi S c Vc(Rdisc), where -Rdisc is the scale length 
of the disc, 1 which can be determined under the assumption 
that the collapsing baryons conserve angular momentum. If 
we ignore the contribution of the halo to the circular velocity, 
for the three cases we find: 



f 5/2 



disc. Mestel 



sc, rigid 



32 ^spin -Rvir 

97T 2 GMl sc f d (Q b /Q M ) 



40A s 2 pln fl vir 



Tdia 



where 



GM 2 m . 



Aspin Rvir f E 



fn 



3f d (n b /n M ) + 5.2fj(n b /n M ) 



(2) 
(3) 
(4) 

(5) 



(see Mo et al. 1998 for the exact expression). At fixed A s 



and f d , Tki, 



sc. rigid 



> T dis 



The total kinetic 



(T) and potential (|W|) energies of the systems, including 
the contribution and stabilizing effect of the NFW halo, in- 
crease with A S pi n due to the halo contribution within i?disc, 
which increases with A S pm. The ratio T/|W|, nevertheless, 
decreases, due to an increasingly dominant halo contribu- 
tion. 

The stability results are summarized in Fig. Q in which 
the maximum spin parameter As P i n ,max for which a disc is 
unstable is shown as a function of the fraction of baryons 



1 Note that both the Mestel disc and the rigid disc are defined 
only for R < R disc (Mestel 1963). 




Figure 1. Bottom panel: maximum spin parameter, A sp in,max! 
for disc instability as a function of the gas fraction ending up in 
the disc. Discs are stable for A sp i n > A sp in,max. Solid line: Mestel 
disc, dashed line: rigid disc, dot-dashed line: exponential disc. 
Upper panel: fraction of unstable discs, for each of the three disc 
models. 



forming the disc, i.e., for every f d , discs are stable for A sp m > 

Aspin, max- 

The distribution of spin parameters found in numerical 
simulations is well fit by a lognormal distribution in A sp m, 
with mean A st 



0.05 and standard deviation ax = 0.5: 



p(X) dX 



2-ira), 



■ exp 



In 2 (A/A) 



<lo\ 



A ' 



(6) 



This function is a good fit to the N-body results of Warren 
et al (1992). Similar results were found in later investiga- 
tions (e.g., Cole & Lacey 1996; Bullock et al. 2001; van den 
Bosch et al. 2002). With this assumption we can estimate 
the fraction of discs subject to dynamical instability, as a 
function of fd, for each of the three disc models (Fig. 0. 



3 STRUCTURE AND EVOLUTION OF 
COLLAPSING GAS 

The unstable conditions described in § 2 are expected to lead 
to runaway infall, provided that the gas remains cooler than 
the local virial temperature as it collapses. For the densities 
expected in pregalactic haloes, the cooling time to < 10 4 
K is much smaller than the dynamical time; this ordering is 
preserved as the collapse proceeds. At the initial disk radius, 
the gravitational potential due to the gas is already appre- 
ciable compared to that of the dark matter. As we will see 
below, collapse leads to a mean gas density profile at least 
as steep as r~ 2 , implying a virial temperature that remains 
roughly constant or increases with decreasing r, whereas the 
dark matter density is expected to increase only as oc r _1 in 
the inner parts of the halo (NFW). The relative dominance 
of gas self-gravity over the dark matter potential thus in- 
creases as the gas collects toward the center, implying that 
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the conditions for large-scale gravitational instability inten- 
sify with decreasing radius. Conditions are therefore ideal 
for the "bars within bars" instability. We will henceforth 
neglect the dark matter. 

The collapse of a self-gravitating, isothermal gas cloud 
has been analyzed in both the nonrotating (Larson 1969; 
Penston 1969; Shu 1977) and rotating but inviscid (Saigo & 
Hanawa 1998; and references therein) limits. In all cases, the 
outer part of the flow evolves toward the density profile of a 
singular isothermal sphere, p oc r~ 2 . Mineshige & Umemura 
(1997), analyzing the case of a shrinking, self-gravitating 
accretion disk subject to an a— viscosity, found analogous 
behavior, i.e., a surface density distribution £ oc r" 1 in the 
outer regions. Using SPH simulations, Englmaier & Shlos- 
man (2004) show that a rotating, self-gravitating gas cloud 
can decouple dynamically from a larger stellar bar, thus ver- 
ifying the basic bars-within-bars picture. They find that the 
shrinking gas bar develops a strong density gradient with 
radius, and an angular pattern speed oc a -1 , where a is 
the semimajor axis of the gaseous bar. These features are 
also consistent with an r~ 2 radial density distribution. 

The p oc r~ 2 behavior can be understood as follows. 
Suppose the initial mass and radius of the cloud are Mo 
and ro, respectively, corresponding to a virial speed vo = 
(GMo/ro) 1 '" 2 ~ lOwio kms -1 . At large radii and early times, 
the evolutionary timescale of the flow is set by the free fall 
time at ro, to ~ ro/vo- At smaller radii the free fall time 
is shorter, but matter is being fed into these radii on the 
much longer timescale to- Therefore, the mass flux through 
these regions must be roughly constant and independent of 
radius, 

M(r) = Mo ~ ^ ~ ^ = °- 2 «io M e yr" 1 . (7) 

to Ct 

On average, the self-gravitating gas at every radius gets rid 
of its angular momentum on the local free fall time, ts(r) ~ 
[GM(r)/r 3 ]- 1/2 , where M(r) is the mass contained within a 
radius r. The mass flow rate is therefore Mo ~ M(r)/tg(r) oc 
[M(r)/r] 3 ^ 2 . Since Mo is constant, we have M(r) oc r and 
p oc r~ 2 . 

This behavior can change if the collapsing gas fragments 
and forms stars. A self-gravitating cluster of collisionless par- 
ticles would increase its velocity dispersion in response to the 
bar potential, thus quenching the instability. Although the 
cooling timescale of the gas is shorter than the dynamical 
time, and therefore violates the Gammie (2001) criterion for 
avoiding fragmentation, we do not expect fragmentation to 
deplete a large fraction of the inflow. This is because the 
circular speed in the potential of the self-gravitating gas re- 
mains roughly constant at r < ro, and corresponds to a 
temperature that is close to the virial temperature of the 
halo. Since we consider only halos with T v i r > 10 4 K, the 
gas never cools very far below T vlI and therefore does not 
form a very thin sheet (i.e., the Toomre [1964] parameter 
Q does not become extremely small). The Jeans mass un- 
der such conditions is only a few times smaller than M(r), 
suggesting that fragmentation will be inefficient. A corol- 
lary of this argument is that the efficient collection of gas 
in the center of a halo might occur only for a relatively 
narrow range of T v i r > 10 4 K, and only under metal- free 
conditions where the formation of H2 is inhibited (these ar- 
guments were discussed at length by Bromm & Loeb 2003) , 



or for gas enriched below the critical metallicity threshold 
for fragmentation (Santoro & Shull 2006). 

As noted in many earlier works, the r~ 2 density pro- 
file cannot persist all the way to the center of the collaps- 
ing cloud (Shu 1977; Mineshige & Umemura 1997; Saigo & 
Hanawa 1998). After a time t, gas collecting at a rate Mo 
will dominate the potential out to a radius ri(t) ~ (t/to)ro- 
An r~ 2 density distribution, with M(r) oc r, is not globally 
self-gravitating inside n , and therefore cannot drive the bars 
within bars instability. In nonrotating collapse (Shu 1997) 
the accumulated gas behaves like a point mass and the den- 
sity at r < n scales as p oc r~ 3 ^ 2 , as in Bondi (1952) accre- 
tion. In inviscid collapse of a rotating fluid (Saigo & Hanawa 
1998), the centrifugal barrier forces the density to be roughly 
constant at small radii. 

Our model involves effective transport of angular mo- 
mentum, and in this respect is closer to the case considered 
by Mineshige & Umemura (1997). Assuming an a— viscosity 
with fixed a, they find that the surface density profile in the 
inner region steepens to E oc r~ 5 / 3 , corresponding to a steep 
mean density distribution p oc r~ S//3 . The corresponding in- 
flow rate scales as M oc r 1 ^ 3 , and the inflow speed, v oc r, 
falls far below the free-fall speed (oc r -1 ^ 3 ). We also expect 
self-similar settling of gas at a mean speed v(r) ~ r/t <§C vo 
for r <JC r\. If we assume that angular momentum trans- 
fer is governed by global self-gravitational instabilities, in- 
stead of an a— viscosity, then the gas must adopt a configu- 
ration where these instabilities are nearly quenched in order 
to transfer angular momentum on a timescale much longer 
than the dynamical time. 

The conditions for global instability depend on the de- 
tails of the gravitational potential as well as the radial dis- 
tributions of density, pressure, and angular momentum. Be- 
cause of the rapid cooling, gas pressure is negligible and we 
expect the system to remain strongly unstable at all radii 
where it is substantially self-gravitating (Shlosman et al. 
1989). The only way to suppress the instability, it seems, is 
for the density distribution to become sufficiently centrally 
concentrated that a large fraction of the gravitational po- 
tential at each r is generated by the gas at much smaller 
radii. This requires the mean surface density distribution 
to steepen to S oc r~ , corresponding to a mean density 
p oc r~ 3 , so that there are roughly equal amounts of mass 
per decade of radius. Note that this is slightly steeper than 
the density distribution obtained by Mineshige & Umemura 
(1997). 

The above argument assumes that angular momentum 
continues to be transported by global gravitational instabil- 
ities, and that fragmentation continues to be unimportant. 
The latter assumption is much less secure at r < n than 
it is further out. If the gas remains isothermal at ~ 10 4 K, 
the Toomre Q— parameter would decrease oc r 1 * 12 and the 
inflowing gas would form a thin disk. Thus, fragmentation 
could seriously hamper the bars within bars instability at 
r < n. If fragmentation is highly efficient, then the inflow- 
ing gas might simply lay down an isothermal stellar potential 
with a constant velocity dispersion ~ vo- Little gas would 
be deposited inside ri(t). 

It seems unlikely, however, that the outcome is this ex- 
treme. Fragmentation should not deplete the gas density 
much below the threshold for local gravitational instability, 
Q ~ 1. Gas with the corresponding density, and a sound 



© 0000 RAS, MNRAS 000, 000-000 



Formation of Black Holes by Direct Collapse 5 



speed ~ vo, would continue to accrete at a rate aMo, where 
a is a viscosity parameter (Shlosman et al. 1990). Even 
where fragmentation is suppressed, local nonaxisymmetric 
gravitational instabilities could continue to drive angular 
momentum transport with an effective a ~ 0(1). Other 
sources of angular momentum transport probably operate 
as well, such as turbulence stirred up by the fragmentation 
and star formation itself. Therefore, although we are not 
able to predict exactly how much gas makes it all the way 
to the central region of the halo, we are probably justified 
in parameterizing the surviving mass flux as aMo, with a 
assumed to be ~ 0(1). 



4 CREATION AND GROWTH OF A 
"QUASISTAR" 

Disk accretion persists as long as the infalling gas is able to 
radiate away the liberated binding energy. Given an accu- 
mulated mass of M,(t) ~ aMot at r < ri, we find that the 
luminosity generated outside a radius r is given by 



L(r, t) ~ aM 



GM* 



2«0 



n(t) 



(8) 



Within a certain radius this radiation is trapped, the pres- 
sure builds up and the gas inflates into a pressure-supported 
cloud, which we dub a "quasistar" . Since rotational support 
no longer dominates, we assume that the self-gravitational 
instabilities are finally quenched. The condition for radiation 
trapping is given by 



L(r,t) >L Edd (t) + 



(9) 



where I/Bdd(i) ~ 4iracGMot/n is the Eddington limit, given 
the appropriate opacity k, for the accumulated mass M*(t). 
p gas and p ra d are the gas and radiation pressure, respec- 
tively, in the quasistar. Once the quasistar mass exceeds a 
few solar masses (i.e., very early in its growth, since the 
mass is growing at ~ 0.2a« 3 Mq yr _1 ), the mean LTE 
radiation pressure exceeds the gas pressure; we will hence- 
forth assume p rad ^> p gas - The interior of the quasistar is 
hot enough to ionize hydrogen, allowing us to assume that 
the opacity is dominated by electron scattering. Substitut- 
ing ri(t) ~ (t/to)ro, we find that the radius of the quasistar 
is time-independent, and is given by 



4tvGc 



1.6 x 10 a«io cm. 



(10) 



Thus, the quasistar that is going to give rise to a supermas- 
sive black hole has a radius of 1 AU, for avf ~ 1- 

Conditions in the interior of the quasistar are ex- 
tremely sensitive to the mass inflow rate, which we ex- 
press through its dependence on a and vo- Denoting the 
quasistellar mass by M* = m* Mq oc t, we find the 
mean density p* ~ 10~ 7 a~ 3 v^ 9 m t g cm -3 , mean pressure 
p* ~ W 6 a^ 4 Vj 12 m 2 erg cm" 3 , and mean temperature (in 
LTE) T» ~ 1.5 x W 5 a' 1 v^rnj 2 K. These estimates justify 
our assumptions that p rad 3> p gas for m„ > a few, and that 
the opacity is dominated by electron scattering. 



5 INTERIOR STRUCTURE AND EVOLUTION 
OF THE QUASISTAR 

The characteristic specific entropy of the matter joining the 
quasistar increases with time, s* = p*/pi^ 3 oc M 2 ^ 3 oc 
i 2//3 . Since hydrostatic equilibrium demands p oc p 2 r 2 oc 



p ' M(r) ' in the quasistellar interior, where M{r) is the 
mass contained within a radius r, we conclude that each 
layer of matter added to the quasistar approximately con- 
serves its specific entropy as the quasistar grows. The posi- 
tive entropy gradient, s(r) oc M{r) 2 ^ 3 oc r 2 ^ 3 , stabilizes the 
quasistar against convection, which would otherwise tend to 
homogenize the entropy. This implies that each layer of the 
stellar interior is compressed by overlying material until the 
radiative diffusion time across the layer, £dift( r ) ~ pur 2 /c, 
is roughly equal to the age of the quasistar. 

The interior structure of the quasistar is therefore 
characterized by a density profile p ~ p„(r/r„)~ 2 , pres- 
sure profile p ~ p*(r/r»)~ 2 , and temperature profile T ~ 
T„(r/r»)- 1/2 . 

These scalings apply as long as radiation pressure ex- 
ceeds gas pressure. However, the decreasing specific entropy 
toward the center implies that the ratio of radiation pressure 
to gas pressure decreases with decreasing r, 



Pgas \ T* / 



1/2 



(11) 



implying that p gas ~ p rad at small enough radii. The gas 
pressure-dominated core comprises the earliest material laid 
down during the growth of the quasistar. It has a radius 
r c ~ r.m^ 1 ~ 1.5 x 10 13 Q« 3 m7 1 cm, temperature T c ~ 
1.5 x 10 5 a _1 u^ 3 m» K, and density p c ~ 10 _7 a~ 3 «^ 9 m 3 g 
cm~ 3 . The core mass is independent of M», a, and «io, and 
is roughly 1 Mq. 

Nuclear burning commences when the core temper- 
ature reaches T c ~ 10 6 — 10 7 K, for a quasistar mass 
ra* ~ (10 — 100)awio. At this point the core density is not 
that dissimilar to densities inside main sequence stars, so the 
burning timescales are likely to be similar as well. The evolu- 
tion timescale due to infall is much shorter, t cv — M*/Mo ~ 
5a -1 uj~ 3 ra» yr, so we do not expect nuclear burning to 
progress very far until the core temperature approaches 
~ 10 8 K, for m* ~ 10 3 q« 3 . At this point the gravitational 
binding energy of the quasistar is ~ 10 13 a _1 w 1 " m* erg g _1 
while the available nuclear binding energy is < 6 x 10 18 m7 1 
erg g _1 , if burning is confined to the core. In order for the 
nuclear energy release to overpower the gravitational bind- 
ing energy of the quasistar as a whole, the mass must satisfy 



m* < 700a 



1/2^3/2 



10 



. Thus, by the time that nuclear reactions 



are able to run to completion, the available energy is prob- 
ably insufficient to seriously affect the outer layers of the 
quasistar. 

It is the ultra-high infall rate, squeezing the core and 
raising its temperature beyond the thermostatic set points 
of thermonuclear reactions, that distinguishes the evolution 
of the quasistar from that of a normal Pop III star. At best, 
nuclear burning can provide a brief hiatus in the contraction 
of the core, which ultimately reaches temperatures > 10 9 K 
where neutrino losses become important. 
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6 CORE COLLAPSE AND INITIAL GROWTH 
OF BLACK HOLE 

Continued compression by infalling matter prevents the core 
from losing energy radiatively and collapsing or becoming 
degenerate. At sufficiently high temperatures, however, neu- 
trino losses lead to core collapse and the formation of a black 
hole. At T c < 10 9 K, the dominant neutrino loss mecha- 
nism is the URCA process, which is ~ 300 times faster than 
photo-neutrino production (Qian &: Woosley 1996; Dutta et 
al. 2004; Itoh et al. 1989; see also Koers & Wijers 2005 for 
a summary of principal rates) . At higher temperatures, pair 
annihilation becomes competitive with the URCA process, 
but because p oc T 3 in the core, both mechanisms scale 
similarly with T c . Therefore, we may approximate the core 
cooling rate by Q c ~ 3 x 10 15 (T C /10 9 K) 9 erg s" 1 cm" 3 , and 
the cooling timescale by t coo \ ~ 4p c /Q c . Setting this equal 
to £ev, we find that the core collapses when m* ~ 3600at>f 
and T c ~ 5 x 10 8 K. 

Details of the collapse depend on the angular momen- 
tum in the core as well as the precise core mass at the time of 
collapse. If angular momentum is initially unimportant, the 
core should collapse at roughly constant temperature. As 
the specific entropy decreases due to cooling, gas pressure 
begins to exceed radiation pressure and neutrino losses are 
dominated by the URCA process. Because the core mass 
is rather low, collapse could get hung up by electron de- 
generacy pressure, but infalling matter from the quasistar 
envelope — which continues to cool via neutrino losses — 
would quickly drive the mass over the Chandrasekhar limit. 
Continued infall similarly circumvents neutron degeneracy, 
with the result that a black hole of a several solar masses 
forms in a few times the core free fall time. 

If the angular momentum of the core and surrounding 
material is too large to permit direct collapse to a black hole, 
then neutrino cooling will lead to a rotationally supported 
disk. As material joins the disk, self-gravity will trigger a 
new round of "bars within bars" instability, which will gen- 
erate additional entropy (enhancing neutrino cooling) as well 
as facilitating collapse by removing angular momentum. 

The amount of matter that falls promptly into the black 
hole depends on the distribution of angular momentum in 
the p oc r~ 2 envelope of the quasistar. The black hole will im- 
mediately swallow all the matter in the quasistellar envelope 
with a specific angular momentum j — Q.r 2 < GMbh/c. At 
one extreme, the specific angular momentum at each radius 
could be a fixed fraction of the Keplerian value, implying 
j(M) oc M as a function of the enclosed mass M(r). In this 
case, the black hole does not grow immediately much beyond 
its initial mass. At the opposite extreme, the internal redis- 
tribution of angular momentum within the quasistar could 
have led to solid-body rotation, implying that j oc M 2 . In 
the latter case, the black hole mass could quickly swallow a 
fraction ~ (GM t /r t c 2 ) x ^ 2 of the envelope's mass, assuming 
that the rotation rate reaches approximately the Keplerian 
value at r*. However, this amounts to only about 20avf M & ; 
therefore the prompt black hole mass is unlikely to exceed 
several tens of solar masses. 



7 SUBSEQUENT BLACK-HOLE GROWTH 

After the initial collapse and prompt accretion phase, the 
growth of the black hole is regulated by angular momentum 
transport. The envelope continues to accumulate mass from 
infall, and the binding energy per unit mass increases with 
time, v*(t) 2 ~ GM»(t)/r,, where r* initially remains con- 
stant. The black hole's gravitational sphere of influence has 
a radius r B H ~ GM B n(t)/v,(t) 2 ~ (M B h/M»)tv If angular 
momentum were unimportant, then the black hole would 
grow at the Bondi rate, M Bondi ~ «* /G, and would swallow 
the quasistar in a free-fall time. Thereafter it would grow at 
the infall rate, aMo. 

Since angular momentum is important, a fraction of the 
binding energy released close to the black hole, eMbhc 2 , 
where e ~ 0.1 is the accretion efficiency, is transported out- 
ward by the torque. If it is not radiated away from close 
to the hole — by neutrino losses, since there is no free sur- 
face from which a wind can emerge — this energy must 
react back on the inflow, slowing down the accretion. Let 
us assume that neutrino losses are negligible. Initially, the 
region affected by the feedback energy is confined to the 
interior of the quasistar. The total energy liberated by the 
time the black hole reaches mass Mbh, £bh ~ eMbhc 2 , 
affects the density profile inside the quasistar out to a ra- 
dius r a ~ E-BH/{p*r 2 v 2 ). Note that the liberated energy is 
trapped inside the quasistar, rather than flowing through 
it in a quasi-steady state, because the radiative diffusion 
timescale at every radius is comparable to the age of the 
quasistar and the growth time of the black hole (§ 5). 

When r a reaches r„, the liberated energy equals the 
binding energy of the quasistar and the latter begins to ex- 
pand. This happens very early in the angular momentum- 
dominated growth phase, when the black hole mass has in- 
creased by an amount AM B h ~ GM 2 /ec 2 r, < O(Mbh). 
To show this, we estimate the rate at which the black hole 
swallows matter from the quasistar envelope. The rate at 
which mass is supplied to the black hole's sphere of influ- 
ence can be taken to be proportional to the Bondi rate, 
M sup ~ 47raBHp(rBH)f»(£)r B H, where p(r B u) is the density 
evaluated at the black hole radius of influence and «bh < 1 
is a parameter that describes the inefficiency of mass cap- 
ture, e.g., due to a finite rate of angular momentum trans- 
port (qbh need not be the same as a). We expect the pres- 
sure and density distributions to be rather flat between tbh 
and r a because of the extra energy injection, and therefore 
use v* to estimate both tbh and the free-fall speed at the ra- 
dius of influence. We also take p(t*bh) ~ p( r a) ~ p» (r,/r a ) 2 . 
Note, however, that even if the density increases oc r~° at 
r < r a , our prescription gives a lower limit to M sup provided 
that (3 < 14/5. 

The rate at which matter actually reaches the black hole 
is suppressed by a further factor, due to the back reaction 
of the energy flux inside the radius of influence (Gruzinov 
1998; Blandford & Begelman 1999; Narayan, Igumenshchev 
& Abramowicz 2000; Quataert & Gruzinov 2000). In the ab- 
sence of a wind that removes energy and/or angular momen- 
tum, the accretion rate is reduced to Mbh ~ (v,/c) 2 M Bup 
(Blandford & Begelman 1999, 2004). We then find that the 
accretion rate is 
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Since v t oc M, 1 ", Mbh oc M* 2 is a steeply increasing func- 
tion of M*. Comparing it to the inflow rate onto the qua- 
sistar, we find that the feedback energy equals the binding 
energy of the quasistar before the black hole mass has dou- 
bled. Thus Bondi accretion, even modified by feedback and 
a finite rate of angular momentum transport, should quickly 
bring the quasistar to the point where its evolution is driven 
by feedback from the black hole. 

The feedback flux does not blow apart the quasis- 
tar, since this would stop the growth of the black hole 
and therefore the feedback. Instead the quasistar expands 
gradually, allowing the black hole accretion rate to ad- 
just so that the feedback energy flux equals the Edding- 
ton limit for the instantaneous quasistar mass, Mbh ~ 
2 x 10~ 3 (e/0.1) _1 (m»/10 5 )M Q yr _1 . The feedback energy 
flux exceeds the Eddington limit for the black hole by a 
factor Af*/A/BH; thus, the black hole grows at a super- 
Eddington rate as long as M» > A/bh- This configuration 
requires most of the feedback flux to be carried by connec- 
tive motions inside the quasistar, since the enclosed mass 
at tbh < r < r* is a steeply increasing function of radius. 
However, one can show that the required convective veloc- 
ity, while larger than the mean inflow speed, is much smaller 
than the local free fall speed at all r. If the quasistar mass 
continues to increase at the constant rate aMo, then the 
black hole mass evolves according to 

M BH (i)~4xl0 5 Q 4(^) _1 (lO^) M ®' (13) 

i.e., M BH oc M 2 . 

To determine the evolution of the quasistar's structure 
in response to feedback, we estimate Mbh using the modified 
Bondi rate discussed above. If we assume the density to be 
roughly uniform within the quasistar (outside run), we have 

Equating this to the Eddington-limited rate (assuming elec- 
tron scattering opacity) and using eq. (I13p and the assumed 
infall rate onto the quasistar, we obtain 

(s)"W(sr-<-> 

Neutrino losses are rapidly quenched by the expan- 
sion of the quasistar. Radiation pressure dominates through- 
out the envelope, and the temperature (in LTE) de- 
creases linearly with time (and with M»), T» ~ 4 x 
10 5 (a BH /0.01)- 1/2 OTj ! o(e/0.1)(m»/10 5 )- 1 K. Within the 
black hole's radius of influence, the pressure varies oc r~ 3 ^ 2 
(not oc r _5,/2 , as in ordinary Bondi accretion, because of the 
feedback), and T oc r -3 / 8 can exceed T„ by a factor as large 
as ( £C > 2 ) ~ 40(a BH /0.01) 3 / 16 ( at ; 1 o)- 3/8 (m t /10 5 ) 3 / 16 , 
close to the black hole. Such temperatures are inadequate 
to produce a significant neutrino flux when the black hole 
grows much beyond its initial mass. 

The effective temperature of the quasistar's photo- 
sphere is also expected to decrease, implying that the qua- 
sistar is unlikely to be a significant source of hard UV ra- 
diation when it has grown beyond ~ 10 4 A/q. Up to this 
point, the rate of production of ionizing photons is very 
high, of order 10 55 photons s _1 ; but since this phase lasts for 
< 10 yr, the total output falls far short of the requirement 



for reionizing the Universe. Similarly, the quasistar produces 
~ 10 50 photons s _1 in the Lyman- Werner band, but can keep 
the molecular hydrogen in its surroundings photodissociated 
only for < 10 5 yr. These estimates correspond to a spherical 
photosphere at r*, but we note that photospheric temper- 
atures could be even lower if the photosphere is strongly 
flattened by rotation. 

The above estimates are valid only as long as T* > 10 4 
K, corresponding to 

M. < 4 x 10 6 a« 3 (^L) M& (16) 

and 

A/bh < 9 x 10 5 ^ 3 (jgy)^ (Jj) M e . (17) 

For the fiducial parameters, the black hole mass at this stage 
is almost as large as that of the quasistar; further growth can 
occur at the Eddington limit of the black hole. However, we 
emphasize the uncertainty in parameters such as qbh and 
a. (We use different fiducial estimates of a and qbh because 
the latter represents viscous transport of angular momentum 
while the former represents a removal of gas from the inflow 
due to fragmentation and star formation.) 

At lower temperatures, the Planck mean opacity (which 
is relevant for calculating the radiation force in LTE, and 
therefore the Eddington limit) becomes very sensitive to 
temperature (Mayer & Duschl 2005), increasing sharply at 
T„ < 10 4 K and then decreasing to several orders of mag- 
nitude below the electron scattering opacity as the temper- 
ature declines further. The sharp decrease in opacity would 
affect the photosphere at an even earlier stage in the qua- 
sistar's evolution. 

Finally, we note the existence of an alternative evolu- 
tionary scenario in which «bh is so small that feedback does 
not regulate the structure of the quasistar. If q?bh is es- 
sentially zero (in practice, <C 10 -6 when m* ~ 10 5 ), then 
the centrifugal barrier forms a wall within the quasistar at 
r ~ vbh- t- m is once again constant and the temperature at 
tbh scales as M*/Mg(j 2 . If this ratio increases with time, fol- 
lowing the initial collapse and prompt accretion phase, then 
neutrino losses remain important within the black hole's 
sphere of influence. We are then justified in assuming that 
self-gravitational instabilities transport angular momentum 
effectively. Provided that nearly all of the liberated bind- 
ing energy is carried away by neutrinos, we deduce that the 
black hole grows at such a rate that Mbh oc M 2 oc t 2 . Adopt- 
ing T(rBH) = 10 9 Tg K as the threshold for rapid neutrino 
cooling, we find that the black hole mass (in solar units) 
grows according to 

m BH ~ 225a-Vo 6 T 9 - 2 (^f ~ 9xl0 4 T 9 " 2 (^±_j (18) 

This rate is not much smaller than the Eddington-limited 
rate, eq. (1131 . in the presence of feedback. Given the con- 
vergence of these two extreme estimates, we are reason- 
ably confident that rapid black hole growth up to masses 
~ 10 4 — 10 6 Mq is possible under the conditions postulated 
here. 

The discussion in §§ 3-7 can be generalized to halos 
with T v i r < 10 4 K. If molecular hydrogen cools the gas down 
to ~ 200 K, then runaway collapse without fragmentation 
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could occur in halos with correspondingly low virial tem- 
peratures. The characteristic infall speed is then wio < 0.2, 
implying inflow rates of a few thousandths of a solar mass 
per year. Nevertheless, a quasistar with a hot dense core 
will eventually develop, and will ultimately collapse to form 
a black hole due to runaway neutrino cooling. The mass of 
the prompt black hole is insensitive to via, and therefore is 
still < 20Mq. However, the mass of the quasistar at this 
stage scales as v\ Q ; thus the quasistar mass is only a few 
times that of the black hole. Moreover, the black hole could 
not reach more than a few thousand solar masses before 
the growth rate becomes sub-Eddington, and limited by the 
inflow of gas into the quasistar. 



8 EVOLUTION OF THE BLACK HOLE 
POPULATION 

How large a population of black holes is likely to result from 
gravitational instability of gas discs in high redshift halos? 
Given the threshold of T v i r > 10 4 K for efficient cooling, and 
therefore for "bars within bars" instability, we can trace the 
co-evolution of the black holes and their hosts. A T v \ r ~ 10 4 
K halo has a mass between 1O 7 M0 and 1O 9 M0 at redshift 
6 < z < 20. The black hole forming in such a host could 
grow at the super-Eddington rate given by eq. 11311 until 
it reaches ~ 1O 6 M0, at which point its mass would ap- 
proach that of the quasistar. Thereafter it could grow, by 
Eddington-limited accretion (or the infall rate, if smaller), 
to an even larger mass. However, the growth of the black 
hole can also be terminated earlier by lack of fuel, i.e., by 
using up all the available gas in the unstable disc (if fd is 
much smaller than unity) and/or if star formation depletes 
the inflow. Therefore, we do not assume that black holes 
grow rapidly to the maximum allowed mass. Indeed, we will 
see below that our mechanism can provide the seeds for all 
present-day SMBHs, even if its efficiency is quite low. 

More massive halos, with T vir > 10 4 K, are probably 
prone to fragmentation and star formation, which would in- 
hibit instability and therefore the formation of a black hole 
by this process. Using the Press-Schechter formalism (Sheth 
& Tormen 1999) , we estimate that the number density of ha- 
los with virial temperature T v i r > 2 x 10 4 K (T v i r > 5 x 10 4 
K) is about 10% (1%) of the density of halos with 10 4 K 
< T vir < 2 x 10 4 K (10 4 K < Tvi r < 5 x 10 4 K). More 
massive halos make an even smaller contribution. Since the 
contribution of halos with T V1I 3> 10 4 K is negligible, we can 
estimate the black hole density by integrating over all host 
halos with T vir > 10 4 K. 

Among all halos with T v i T > 10 4 K, only those with a 
low enough spin parameter (given fd) will host a disc unsta- 
ble to "bars within bars" instability. Assuming a seed black 
hole mass of 2OM0, we plot the comoving mass density of 
black hole seeds in Fig. |2] The mass density is small, but 
this process can nevertheless seed most of the systems which 
will evolve into the local galaxies where SMBHs have been 
found. For a given local halo, the Extended Press-Schechter 
formalism can be used to estimate the average number of 
progenitors with T v i r ~ 10 4 K at z > 10. The probability 
of black hole formation depends also on the amount of gas 
which condenses to form a disc (see § 2). If the fraction of 
gas typically ending up in the disc is fd ~ 0.5, we find that 



this process needs to operate only until z ~ 18 in order 
to supply seeds to all present-day halos of Milky- Way, or 
larger, size. If black hole formation proceeded to z ~ 14, 
then all halos with mass > 10 M@ today could have been 
seeded (see Barth, Greene & Ho 2005). If fd was lower, the 
black hole formation process would have had to continue for 
longer, in order to form seeds in the progenitors of most lo- 
cal galaxies. These constraints would be eased by an early 
generation of black hole seeds, formed in a small fraction of 
halos with XVit > 200 K, before H2 is photodissociated by 
the first stars. 

The process of seed formation can be widespread 
enough to account for subsequent SMBH evolution, even 
if the subsequent black hole growth (§ 7) is inefficient in 
most high-redshift halos. We note, in fact, that black holes 
with masses well below ~ 1O 6 M0 are expected in local faint 
AGN (Barth, Greene & Ho 2005). The growth of black hole 
seeds up to ~ 10 4 — 1O 6 M0 cannot therefore happen in all 
high-redshift systems. 

We can estimate the upper limit to the total black hole 
mass density predicted by our model by assuming continu- 
ous formation of seed black holes, and adopting eq. 1131 1 to 
estimate their early growth. We let black holes grow until 
either they have consumed all the gas, or the infall of mat- 
ter onto the black hole greatly decreases, at Mbh ~ 10 6 Mq. 
The comoving density of halos is estimated using the Sheth 
& Tormen (1999) halo mass function. The black hole mass 
density shown in Fig.|2]must be compared to the density that 
we observe at low z, i.e., pbh(z = 0) w 3 — 5x 10 5 M Q Mpc" 3 
(Yu & Tremaine 2002; Aller & Richstone 2002; Marconi et 
al. 2004; Fabian & Iwasawa 1999; Elvis, Risaliti & Zamorani 
2002) and p B n{z = 3) « 4 - 5 x 10 4 M Q Mpc" 3 (Merloni 
2004). At z = 6, a lower limit to the black hole density 
is obtained by integrating the observed luminosity func- 
tion, pbh(z = 6) « few x M Mpc" 3 (Fan et al. 2004). 
This density, however, includes only black holes with masses 
> 10 9 M©. To obtain a rough estimate of the total black hole 
density at z = 6, we adopt Soltan's (1982) argument. We ex- 
trapolate the luminosity function up to z = 6, assuming the 
redshift dependence given by Richards et al. (2005), inte- 
grate the emitted quasar luminosity, and convert it into a 
black hole density by normalizing to the z = SMBH den- 
sity. We do not correct here for obscured quasars, or assume 
a radiative efficiency, and therefore consider the density that 
we find, pbh(z = 6) « 1 — 3 x 1O 3 M Mpc" 3 , more of an 
indication than a robust estimate. 

The densities shown in Fig.|5|must be regarded as upper 
limits to the black hole density, as we have not included any 
effect that can interrupt or disturb black hole formation and 
growth. We discuss in the next section how the hierarchical 
framework for structure formation can modify this simple 
picture. We also recall that, eventually, efficient star forma- 
tion occurs in these halos, competing for the gas supply and 
possibly limiting the mass available for black hole accretion. 



9 DISCUSSION AND CONCLUSIONS 

We have presented a scenario for the accumulation of gas 
in the centers of dark matter halos with T v j r > 10 4 K, the 
initial collapse of the gas to form a seed black hole, and 
the subsequent early growth of a supermassive black hole. 
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Figure 2. Lower panel: comoving density of 20Mq black hole 
seeds as a function of redshift. Upper panel: comoving density of 
black holes, assuming continuous formation and growth of seeds 
to 10 6 Mq according to eq. 1131 . starting from z = 20. Solid line: 
Mestel disc, dotted line: rigid disc, dot-dashed line: exponential 
disc. The upper set of lines assumes fd = 0.5, the lower set as- 
sumes = 0.1. 



This mechanism can lead naturally to the super-Eddington 
growth of black holes up to masses ~ 1O 6 M0, as early as 
redshifts 10 — 20. Given additional growth to ~ 1O 9 M0 at 
close to the Eddington rate, the model can account for the 
population of quasars observed at z ~ 6 (Fan et al. 2004). 
Even without significant growth after the formation phase, 
this mechanism could produce the seeds for all supermassive 
black holes inferred to exist in the local universe. 

We argue that global self-gravity triggers the "bars 
within bars" instability (Shlosman et al. 1989, 1990), un- 
der certain conditions, as gas forms a rotationally supported 
thick disk in the center of the halo. On scales much smaller 
than the disc radius, and times shorter than the free-fall 
time, quasi-steady inflow is a better representation of the in- 
fall than a monolithic collapse. Local, or quasi-local sources 
of "viscosity", such as those due to magnetic fields, turbu- 
lence, or radiation drag, are not required to transport the 
angular momentum that inhibits black hole formation. In 
metal-free halos with little molecular hydrogen, this behav- 
ior is possible once the virial temperature exceeds ~ 10 4 
K (Oh & Haiman 2002). Under these conditions, gravita- 
tional instabilities can transport angular momentum effec- 
tively from scales of several kiloparsecs down to scales ini- 
tially as small as ~ 1 AU, at a fraction of a solar mass per 
year (for a characteristic infall speed no ~ 10 km s~ ). We 
suggest that inflow is most efficient in halos where T V [ T does 
not exceed a few times 10 4 K, since fragmentation of the 
infalling gas is unlikely to be efficient in this case (Bromm 
& Loeb 2003). As the mass in the center builds up, global 
instabilities may be quenched in the inner regions, but local 
gravitational instabilities could continue to drive a substan- 
tial inflow, even if a certain amount of star formation occurs. 



Instability occurs only where the halo spin parameter, 
A S pi n , falls below a threshold value that depends on fd, the 
fraction of gas that forms the disc. For fd > 0.5 the thresh- 
old is comparable to the mean spin parameter predicted by 
simulations, and > 20% of all T v i r ~ 10 4 K halos should ex- 
hibit instability. Even a value of fd as small as 0.1 leads to 
instability in > 1% of halos and a significant seed population 
of black holes. 

In halos with the low angular momentum required to 
trigger black hole formation, the rapid infall of gas leads 
to a mass accumulation much larger than that expected in 
a mini-halo with an average spin parameter. The forma- 
tion of a "standard" Pop III star (Bromm, Coppi & Larson 
1999, 2002; Abel, Bryan & Norman 2000) is therefore sup- 
pressed in favor of a massive "quasistar" (§ 4). We suggest 
that the much smaller mini-halos (with virial temperature 
well below 10 4 K) that form the first stars are distinct from 
the halos that form the seeds of supermassive black holes, 
although the former may be precursors of the latter. Pho- 
todissociation of molecular hydrogen, possibly by a small 
population of Pop III stars, would suppress fragmentation of 
the infalling gas. It is therefore possible that the formation 
of seed black holes follows an earlier epoch of star forma- 
tion, as the "quasistar" itself is not a significant source of 
photo-dissociating photons for long. The epoch of black hole 
formation must happen early enough, however, that the Uni- 
verse is not highly metal enriched — later episodes of star 
formation would enrich the environment of seed black holes. 

The most important conclusion of our model is that 
the "quasistar" formed by the accumulating gas has a low- 
entropy, gas pressure-dominated core surrounded by a ra- 
diation pressure-dominated envelope. As matter piles onto 
the quasistar, the core is squeezed until its temperature ap- 
proaches 10 9 K (typically when the envelope mass reaches 
a few thousand Mq). Cooling by thermal neutrinos then 
leads to core collapse and the formation of a seed back hole 
of ~ 10 — 20Mq . This is a novel application of neutrino cool- 
ing, which has been invoked previously in connection with 
hyperaccretion onto neutron stars in supernovae (Colgate 
1971; Chevalier 1989; Houck & Chevalier 1991) and com- 
mon envelope binaries (Chevalier 1993; Brown, Lee & Bethe 
2000), and onto black holes in gamma-ray bursts (Narayan, 
Paczynski & Piran 1992; Woosley 1993; Popham, Woosley 
& Fryer 1999; Narayan, Piran & Kumar 2001). It is difficult 
to set up the necessary conditions for efficient neutrino cool- 
ing, since radiation pressure generally prevents the accretion 
rate from reaching the required level from below, unless the 
viscosity parameter a is extremely small (Chevalier 1996). 
Previous models have circumvented this problem by invok- 
ing strong radiation trapping in a steady-state (or slowly 
varying) accretion flow. In our case the inflow sets up favor- 
able conditions in a time-dependent fashion by establishing 
a steep positive entropy gradient in the quasistar, with only 
mild radiation trapping. 

The subsequent evolution of the black hole can be very 
fast, with growth to more than a million solar masses possi- 
ble in less than a Salpeter timescale. Even taking account of 
strong energy feedback driven by angular momentum trans- 
port, we conclude that black holes can accrete at the Ed- 
dington rate for the quasistar mass, which exceeds the Ed- 
dington rate for the black hole by a factor M»/Mbh- For 



© 0000 RAS, MNRAS 000, 000-000 



10 M. C. Begelman, M. Volonteri & M.J. Rees 



steady infall onto the quasistar, this corresponds to a black 
hole mass increasing with time as Mbh oc t 2 . 

If black hole growth (eq. MM! proceeded undisturbed in 
all halos satisfying the instability criterion with T vlv > 10 4 
K, then the total mass density in supermassive black holes 
would become comparable to the local one already at high 
redshift. However, a number of effects can limit this initial 
phase of rapid growth. Limitations intrinsic to the halo in- 
clude the overall mass supply that participates in the infall, 
as well as removal of matter from the inflow by star forma- 
tion. Moreover, the halos and their embedded black holes do 
not grow in isolation. Halos susceptible to the "bars within 
bars" instability represent high peaks in the field of den- 
sity fluctuations (Tegmark et al. 1997; Madau et al. 2004). 
Therefore, they should experience an enhanced number of 
major mergers with respect to "average" halos at the same 
redshift. Halo major mergers can modify our basic results in 
two ways. First, cosmological simulations show that the spin 
parameter of a halo typically increases after a major merger 
(Vitvitska et al. 2002). On the other hand, the spin param- 
eter decreases after a long series of minor mergers. Major 
mergers, therefore, should delay the triggering of instabili- 
ties, at least until a sufficient number of minor mergers has 
lowered the spin parameter again. Second, a major merger 
could destroy the coherence of the "bars within bars" pro- 
cess. By interfering with the infall of matter onto the qua- 
sistar, a violent encounter could hasten the depletion of the 
mass supply well before the upper limits discussed above are 
reached. Such a disturbance is unlikely to modify the inte- 
rior structure of the existing quasistar or suddenly stop the 
growth of the black hole, however, since these involve only 
the very core of the system. 

Seeding of larger halos by black hole mergers could also 
be limited by the "gravitational rocket" effect, the recoil due 
to the net linear momentum carried away by gravitational 
waves in the coalescence of two black holes (Madau et al. 
2004, Haiman 2004, Yoo & Miralda-Escude 2004). The recoil 
velocity still has large uncertainties, but can easily exceed ~ 
100 km s - , comparable to the escape velocity from shallow 
halo potentials. Volonteri & Rees (2006) estimate that up 
to 50% of black holes merging in high-redshift halos can be 
ejected due to the gravitational rocket effect. 

Despite these potential sources of inefficiency, the mech- 
anism we have outlined could be the principal route leading 
to supermassive black hole formation in galactic nuclei. The 
main elements of the model — particularly the cascading in- 
fall via "bars within bars" instability, and the formation and 
evolution of the quasistar with runaway neutrino cooling — 
should be testable via numerical simulations. We hope that 
such simulations can be undertaken shortly. 
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